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Abstract 

We study the Higgs sector of the next-to-minimal supersymmetric standard 
model (NMSSM) with explicit CP violation at the one-loop level, where the ra- 
diative corrections due to the quarks and squarks of the third generation are taken 
into account. We expect that, within a reasonable region of the parameter space 
of the present model, at least one of five neutral Higgs bosons may be produced at 
the future e^e~ International Linear Collider (ILC) with y/s = 500 GeV, with cross 
section larger than 12 fb, 15 fb, and 1.5 fb, respectively, via the Higgs-strahlung 
process, the WW fusion process, and the ZZ fusion process. We find that the 
effect of the CP phase in the present model yields significant influences upon the 
production cross sections of the five neutral Higgs bosons. We also study the decay 
modes of the five neutral Higgs bosons to find that their decay widths are similarly 
affected by the CP phase. Some of the decay modes in the present model behave 
differently from those of the Standard Model. 
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I. INTRODUCTION 



These days, supersymmetry (SUSY) is one of the most important topics for both the- 
oretical and experimental high energy physics. The phenomenology of SUSY has been 
investigated quite exhaustively recent years with respect to the future high energy exper- 
iments at the Large Hadron Collider (LHC) or at the International Linear Collider (ILC). 
From a phenomenological point of view, SUSY must be broken in order to be realized in 
nature [1]. The breaking of SUSY may be accomplished by several methods, one of which 
is by introducing soft SUSY breaking terms. 

The exact form of the soft SUSY breaking terms are not precisely determined. They 
may differ from one model to another. In particular, the coefficients of these SUSY 
breaking terms and their associated terms may either be real or complex. The complex 
phases that may arise from these soft SUSY breaking terms in a supersymmetric extension 
of the standard model (SM) can be the source of CP violation, through the mixing between 
the CP even and the CP odd states, since there are at least two Higgs doublets in the 
model in order to give masses to the up-quark sector and the down-quark sector separately 
[2]. 

It is already known that, at the tree level, neither explicit nor spontaneous CP violation 
is possible in the Higgs sector of the minimal supersymmetric standard model (MSSM), 
which has just two Higgs doublets, because any complex phases in its Higgs sector can 
always be eliminated by rotating the Higgs fields. Even at the one-loop level, spontaneous 
CP violation is disfavored because it requires a very light neutral Higgs boson, which has 
already been ruled out by experiments. On the other hand, explicit CP violation is viable 
in the MSSM at the one-loop level since the radiative corrections due to the loops of 
relevant particles, such as quarks and squarks, yield the mixing between the CP even and 
the CP odd neutral Higgs bosons [3-10]. 

The next-to-minimal supersymmetric standard model (NMSSM), with an additional 
Higgs singlet besides two Higgs doublets in its Higgs sector, may not accommodate spon- 
taneous CP violation at the tree level because of vacuum stability [11]. However, the 
scenario of explicit CP violation is possible for the NMSSM at the tree level. By assum- 
ing the degeneracy of the stop quark masses in the Higgs sector of the NMSSM, large 
explicit CP violation may be realized as the vacuum expectation value (VEV) of the neu- 
tral Higgs singlet approaches to the electroweak scale [12]. The tree-level Higgs potential 
of the NMSSM may have only one physical CP violating phase, either by redefining the 
phases of the three Higgs fields or by applying two CP-odd tadpole minimum conditions 
in the explicit CP violation scenario. 

In previous articles, we have analyzed the phenomenology of the explicit CP violation 
in the Higgs sector of the NMSSM at the 1-loop level, where the radiative corrections due 
to quarks and squarks of the third generation are taken into account [13]. There can be 
several complex phases which allow the explicit CP violation in the Higgs sector of the 
NMSSM at the one-loop level: Two complex phases may arise from the stop and sbottom 
quark masses that are associated with quadratic coupling A and the trilinear SUSY break- 
ing parameters. At and Aj,. Additional complex phases from the chargino and neutralino 
sectors may also contribute to the CP mixing among the scalar and pseudoscalar Higgs 
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bosons. 

The presence of explicit CP violation may affect the Higgs phenomenology such as 
the masses of the neutral Higgs bosons, their productions and decays. In high-energy 
e"'"e~ collisions, the dominant production channels for neutral Higgs bosons are the Higgs- 
strahlung process, the WW fusion process, and the ZZ fusion process. Within the context 
of explicit CP violation, the prospects of discovering the neutral Higgs bosons in high- 
energy e~^e~ collisions have been studied in a general two-Higgs doublet model [14] and 
in the MSSM [15,16]. In particular, in the MSSM with explicit CP violation, where the 
CP mixing among the CP even and CP odd states is maximized, the production cross 
sections for the neutral Higgs bosons in e^e~ collisions at ^/s = 500 and 800 GeV are 
calculated [16]. Meanwhile, the decays of the Higgs bosons in the MSSM with explicit 
CP violation have also been studied in the literature [5,10]. 

In this article, we study the Higgs sector of the NMSSM at the one-loop level in the 
explicit CP violation scenario, where the quarks and squarks of the third generation are 
taken into account. The two CP-odd tadpole minimum conditions are obtained at the 
one-loop level. There appear three physical phases as the possible sources of explicit CP 
violation in the NMSSM. We establish and explore a reasonable region in the parameter 
space of the present model, including the CP phases, in order to examine how the upper 
bound on the masses of the neutral Higgs bosons are dependent on the parameters. 

Further, the production cross sections, and decay rates as well, of the neutral Higgs 
bosons are calculated for the future e~^e~ ILC with ^/s = 500 GeV, within the established 
parameter region, via the Higgs-strahlung process, the WW fusion process, and the ZZ 
fusion process. We study the effects of the CP phases upon the production cross sections 
of these neutral Higgs bosons and their decay rates. We also obtain the lower bound on the 
production cross sections of these neutral Higgs bosons within the established parameter 
region, which allow us to expect that the neutral sector of the NMSSM with explicit CP 
violation may be tested at the future ILC with ^/s = 500 GeV. 

II. THE HIGGS SECTOR 

The Higgs sector of the NMSSM consists of two Higgs doublet fields Hi = 
H2 = {H^^H^), and a neutral Higgs singlet field [17,18]. Their hypercharges are -1/2, 
1/2, and 0, respectively. For the fermion matter fields in this model, only the quarks of 
the third generation are taken into account. Then, the superpotential of the model may 
be expressed as 

W = htQH^f^ + hbQHih]^ + XNH^eH2 - ^N' , (1) 

where hq{q = t, b) are the Yukawa coupling coefficients, Q is the chiral superfield contain- 
ing the left-handed t, b quarks, g|j (g = t, b) are the charge conjugates of the right-handed 
t,b quarks, and e is an antisymmetric 2x2 matrix with 612 = 1. Note that all coupling 
coefficients are dimensionless. 

As is well known, the parameter A is introduced into the NMSSM to avoid the so-called 
/i-problem in the MSSM, since Ax corresponds to the /i parameter in the superpotetial of 
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the MSSM, where x is the VEV of A^. The global U(l) Peccei-Quinn symmetry is explicitly 
broken by the presence of the cubic term proportional to k in the superpotential. Without 
the k-teim, a massless pseudoscalar Higgs boson would appear beacuse the Peccei-Quinn 
symmetry is recovered at the tree level. Thus, the two dimensionless parameters A and k 
are essential in the NMSSM. 

The Higgs potential at the tree level may decompose into three parts as [18] 

Vo = Vd + Vf + Vs , 

where 

Vn = ^fiHlm + HlaH^r + ^fm\' - \H,\')\ 

Vf = |A|2[(|ifi|2 + \H2\'')\N\^ + l^fe^2|'] + lA^HiVl^ - {\k*H^eH2N*^ + H.c.) , 
Vs = ml^\H,\'' + mljH2\^ + ml\N\^-{XAxH^eH2N + ^kAkN' + R.c.) , (2) 

with gi and g2 being the U(l) and SU(2) gauge coupling constants, respectively, and 
a = (cr^, 0"^, cr^) are the Pauli matrices. The soft SUSY breaking Vs has three soft masses 
TTT'Hi, '^H2^ ^-^d mAT, and two additional parameters Ax and A^, both of mass dimension. 

In general. A, /c, A\, and Ak in Vq can be complex. Among them, \A\ and kAk can 
be adjusted to be real and positive by redefining the phases of Hj^eH2 and A^. Thus, at 
the tree level, /ce*'^* and Ae**^^ can in general be complex. Consequently, the tree-level 
Higgs potential may have at most only one physical phase, which may be chosen to be 
= 0A - 0fc in \k*. 

The radiative corrections at the one-loop level to the above tree-level Higgs potential, 
obtained by employing the effective potential method, is given as [19] 



niMt 
647r2 



A^-2 



(3) 



where the subscript / stands for the quarks and squarks of the third generation, rig = —12 
and rig = 6 are respectively the degrees of freedom for quarks and squarks, and A is 
the renormalization scale in the modified minimal subtraction scheme. In the above 
expression, note that quarks enter with a negative sign while squarks enter with a positive 
sign, and M.f depends on the neutral Higgs fields, where the tree-level masses of the 
relevant quark and squarks should be used as input. The radiatively corrected Higgs 
potential at the one- loop level is given hy V = Vq + Vi. 

After the electroweak symmetry breaking, neutral Higgs fields H^, H2, and A^ may 
acquire complex VEVs vi, f2e*^, and xe''^, respectively, where vi is adjusted to be real. 
The tree-level masses of the quarks of the third generation are given as rrit = htv sin (3 



and rrii, = hhVcos(3, where tan/3 = V2/V1 and v = yvf + V2 = 175 GeV. 
The masses of the squarks of the third generation are 



m~_^ t2 ~ + "^t -f 'n^'t^J + A^x^ cot^ (3 + 2XAtX cot (3 cos ift , 

h ~ + "F f^b^/ + ^^^^ t^^^ + 2Ay46a; tan f3 cos 0^ , (4) 
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where mg is the soft SUSY breaking mass, and At and are the trihnear SUSY breaking 
parameters of the mass dimension. Two CP phases 0f and 0^ are given by 0j = ^ + 5 + 
0A + <f>At cLiid (f>b = 9 + S + (f)x + 0^6 5 where (px and ipAgiQ = t, b) come from the complex 
phases of A and Aq{q = t,b), respectively. Note that the masses of squarks of the third 
generation are different from the masses of the quarks of the third generation, due to 
the SUSY breaking terms, which account in general for the mass splitting between an 
ordinary particle and its superpartner depending on the SUSY breaking scale (1-2 TeV). 
For simplicity, we reexpress the Higgs multiplets via a unitary transformation as 



Hi 



vi + 5*1 + z sin (3 A 
sinpC+* 



^ COS/5C+ \ (5) 

^ I {v2 + S2 + icospA)e'' ) ' 
N = ( {x + X + iY)e'^ 



where 5*1, 5*2 and X are the scalar Higgs fields, A and Y are the pseudoscalar Higgs 
fields, and is the charged Higgs field. Note that the phases are explicitly shown in the 
expression. 

As we allow CP violation in the Higgs sector of the NMSSM at the one-loop level, 
there are several complex phases: 0^, 4>\y (pAx^ 'pAk: 'pAt, and (pA^ These phases arise 
respectively from k, A, Ax, A^, At, and A^. In addition to these phases, we have two more 
complex phases, namely, 6 and 6. Thus, we have eight complex phases in the radiatively 
corrected Higgs potential. The two CP-odd tadpole minimum conditions, with respect to 
A and Y, are obtained as 



3/j2 

kxsm(f) + Ax sin 01 + j^^t sin 0t/(m?^, m|) 
3/^2 

— 2kv^\x sin 2(3 sin + AxXv"^ sin 2(3 sin 0i — 2kAkx'^ sin i 
3h't 



' iQjj^2 t ^ -rtj \ ti^ t2' 

+ — ^AAt;2sin2/3sin0,/(m?^,m|) , (6) 



where the dimensionless function f{ml,my) is defined by 



f{ml,ml 



1 



fm?, — 



i2 /m2 



mt 2 1 "^t 



log — - log — - 



+ 1 , (7) 



and = 0A - 0fc - 2(5 + 6, 0i = 0a + 0Aa + d + 6, and 02 = 0fc + <pAk + 35, with 
0f = ^ + (5 + 0A + (pAt and 06 = 6* + 5 + 0a + 0Af Note that here at the one- loop level is 
different from the phase previously defined at the tree level. After imposing these two 
CP-odd tadpole minimum conditions, we are left with only three phases and (pt and 0b 
in the radiatively corrected Higgs potential V at the one-loop level. These three phases 
are the physically free parameters. 
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For the masses of the five neutral Higgs fields, we differentiate V with respect to the 
neutral Higgs fields, in order to obtain the symmetric 5x5 mass matrix M for the neutral 
Higgs bosons. The matrix elements are calculated in the (5*1, 5*2, A, X, y)-basis as 

Mil = Mfi + {niz cos + m\ sin^ j3, 
M22 = MI2 + {mz sin jS^ + m\ cos^ ^, 
M33 = MI3 + m\ , 

M44 = MI4 H -ml sin^ 2/5 Xk sin 2/3 cos </> + (2kxY - kxA^ cos </)2 , 

4x"' 2 

M55 = M55 + —^rnl sin^ 2/3 H Xk sin 2/3 cos + 3kxAk cos 02 , 

4x"' 2 

M12 = Mf2 + (2A\;^ - ^1 - "1^) sin /? cos /3 , 

Mi3 = Mi^3, 

Mi4 = Mi4 m\ sin^ /3 cos /3 — vXkx sin /? cos + 2vX'^x cos /? , 

X 

M15 = Ml 5 H — 3AA;fa;sin/3sin0 , 
M23 = MI3 , 

M24 = M24 mi sin /3 cos^ /3 — vXkx cos /3 cos + 2t>A^x sin /? , 

X 

M25 = M25 — SXkvx cos /3 sin , 
M34 = MI4 + AA;fxsin0 , 

M35 = M35 H — sin P cos /9 — SvXkx cos ^ , 

M45 = MI5 + 2AA;t;^ sin 2/3 sin , (8) 

where m\ is introduced for convenience as 

2 Ax(/1a COS01 + /cxcos0) 3m'fAtXx cos (j)t r-( 2 2\ 
~ sin /3 cos /3 167r2^;2 sin^ /3 cos /3 ^"^'^ ' "^^^ 

SmgAfeAa: cos 0fc 2 2n /q^ 



cos^ /3 sin /3 

generations. Explicitly, the contributions due to top and stop quarks are given as 
SmfX^x^Al g{ml, m?) 



and M?- = Mj- + M^- are the matrix elements due to the contributions of the third 

''J *j 



M- 



87r2f2sin^/3 (m?_^ — m?^)^ ' 

3<^?A| (7(m|, m|) Sm^A^A^^ log(m|/m|; 



STT^f^ sin^ /3 (m?_^ — m?^)2 An'^v'^sm^P (mj — m; 



22 Q^2„.2r.i^2 _2 ^2 /l^2„,2 oiv.2 ^ ^_2 _ ^2 ^ 

^87rVsin2/3 ^1^ ^ ' 

i _ 3m^A Vv4^ sin^ 0t m|) 
~ 87rVsin^/3 (m| -m|)2 ' 
Sm^A^Al ^(m|, m|) 



M; 



87r2tan2/3(m| -m|)2 ' 
, _ 3m^AM2sin2 0i c/(m|, m|; 



87r2tan2/3 (m|-m|)2 



_ 3mtXxAtAi^Ai^ gjml, m|) Sm^AxA,-^ log(m|/m: 



^"^ 87r2t)2 sin^ /5 (m| — m?^^ Svr^ti^ sin^ /3 (m| — m?J 

_ 3m^AVAfA,-^sin0, ^Kv m|) 

87rVsin3/5 (m| -m|)2 ' 



^ 3m^A^xA| ^K, m|) 3m^A^xcot/3 



_ _ 3m^A2xAAiiSin0t (?(m|, m|) 



J _ 3m^XxAfAt^ sin 0j , m| ) Sm^AxA^ sin 0^ log(m| /m?^ ) 



M 



SmnAtAr Ar ^7(m? m?) / 3mnAr log(m?/mn 



87r2f sin /? tan /5 (m? — m? )2 871^1; sin /5 tan /3 (m? — mj ) 



M; 



t2 

, _ 3mfAA2A,-^sin0f ^(m|, m|) 3m^AAsin(/>f log(m|/m|) 

87r2t; sin /? tan /3 (m|_^ — 87r2w sin /3 tan /? (m?^— m|^) 

.2 ^2 \ 



^ ^ _ 3m^A^xAA^-^sin0t ^K,. ^t^' 

87r2?;sin2/5tan/? {m\-m\y' 
^ 3m^A^xA2sin^(/>, ^(m?^, m|) 
3^ 87r2?7sin2/3tan/3 (m| - m|)2 ' 
, ^ 3m^A^AA,-^sin0/^(m|i m|) 

87r2tan2/? - ^ ' 

and the contributions due to bottom and sbottom quarks are given as 

^ 9{ml, ml) ^ 3mtA,A~,^ log(m|/mgJ 

87r2t>2 cos^ /5 (m? — Ati'^v'^ cos"^ (3 (mr — ni^ ) 

/r^ .. log' 



87r2f 2 cos^ /5 I mi 

I2 9i^b "^V 

Stt'^v'^ cos^ /3 [m^ — m^^Y 



,2 ^2 



, ^ 3m^A^xMg sin^ 0fe 
3^ 87rVcos4/? (m? - m? )2 ' 



3mfA^A? ofm? , m? 

j^b _ " 62 ^ ^ fei ' 62 



87r2 cot2 /5 (m? - mn2 ' 
^02 61' 



3m^AM2sin^0b ^K^, 

87r2cot2/5 (m? - m? )2 
^02 01 ' 



^ 3m^AxA,Aj;^A^^ ^(mf^, m?^) 3m^AxA^^ log(»-^iy< 

12 o^2„,2„^o2/Q /'^2 ^2 \2 ~'~ o^2„,2 „^o2 /5 /'^2 ^2 



M- 



87r2t;2 cos2 /5 (m? — m? )2 87r2t;2 cos2 /5 (m? —mi) 

' ^ 02 bi' ' ^ 02 01 ' 

3m^AxA2A^^ sin 06 ^(m?^, 3m^AxAsin0fe log(m?ym?J 



b 

""^^ 87r2t>2 cos^ /3 (m? — )2 8-7r2t>2 cos'^ /3 (m? — ' 
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""^^ Svr^f cos /9 cot /? (m? — m? )2 Svr^f cos /3 cot /3 (m? — m? ) ' 
_ 3rntXAlA~,_^ sin/, ^K|, m|) Sm^AA, sin ^/logK^mf; 



^vr^f cos /9 cot (m? — m? )^ Svr^f cos /9 cot /3 (m? — m? ^ 

2 N ^ 



3m^AVAA^^sin0fe 



STT^ti^ cos^ /5 (rrii — rrii )^ ' 
SmfX'^xAl q(ml , ml) SmlX'^x tan B 
87r\;cos/5cot/5(m?^ -m?^)2 Syr^t^cos/? '^^ "^i' ^2^' 
3m^A2xAA^^sin06^(m?^, m|) 



Svr^t) cos /5 cot /5 (m? — m? )2 

' ' ^02 01' 

^ _ Sm^A^xAA^^sin^ft ^(m?^, m|) 



2^ 87r2t;cos2/3cot/3 (m? - m? )2 

^02 01' 

Sm^A^xA^ sin^ 0^, g{ml , m|) 



Svr^i; cos2 cot (m? — rrij 
^02 01' 



Sm^AMftAi sin 06 



^5 87r2cot2/? (m-2 - m? )2 ' ^^^^ 

' ^02 Ol' 

where 

Ai^ = At cos 0t + Ax cot /3 , 
A^2 = At + Ax cot P cos 0t , 
A^^ = + Axtan/3cos0fe , 

A^^ = y4ftcos0t + Axtan/3 , (12) 
and the dimensionless function g{m'^,my) is defined by 

2 2^ "1+^1 "^^^o nQ^ 
^(m^,m) = ^ + 13 

— rriy 

Note that since these matrix elements are calculated at the one-loop level, they naturally 
contain /(m?^,m?_^) or /(m?^,m?^). In the above expressions, we use the two CP-odd 
tadpole minimum conditions to eliminate most of them. 

The matrix elements that are responsible for the scalar-pseudoscalar mixing are M13, 
Mi5, M23, M25, M34, and M45. The physical neutral Higgs bosons and their squared 
masses are given respectively by the eigenvectors and the eigenvalues, denoted as (z 
= 1-5), of the mass matrix. We sort these five neutral Higgs bosons in the increasing 
order of their masses such that, for example, m\_^ is the smallest eigenvalue and hi is the 
lightest neutral Higgs boson. 

Finally, the charged Higgs boson mass at the tree level is obtained as 

2Ax 

= mL — A^f^ H -(A^cos^i + kxcoscp) . (14) 

sin2p 

where is the W boson mass. Note that, unlike the case of the MSSM, the tree-level 
mass of the charged Higgs boson may be either heavier or lighter than W boson, according 
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to whether the second term is smaller than the third term. Also note that the charged 
Higgs boson mass contains cos 0, even though there might be no CP mixing in the charged 
Higgs sector. 

III. NUMERICAL ANALYSES 

For numerical analyses, we take sin^ 6^ = 0.23 for weak-mixing angle, Gp = 1.166 x 10~^ 
for the Fermi coupling constant, Tz = 2.48 GeV for the total decay width of Z boson, 
mz = 91.187 GeV for the mass of Z boson, rriw = 80.423 GeV for the mass of W boson, 
rrit = 175 GeV for the mass of top quark, and = 4 GeV for the mass of bottom quark. 
In the NMSSM with explicit CP violation, we have as many as thirteen free parameters: 
A, 0, 0f, (j)b, niQ, tan/3. A, k, Ak, x, At, Af,, and A\. Among them, we set A = 300 GeV, 
0i = 06, and At = Af, for simplicity. Then, we are left with ten free parameters 0, (pt, fnq, 
tan/3. A, k, Ak, x. At, and A\. However, we will use rriA in Eq. 9 instead of A\, since it 
is the lighter pseudoscalar Higgs boson mass in CP-conserving limit. 

We explore a reasonable region of the parameter space of the model, whose boundaries 
are set as < 0, 0f < vr, 2 < tan /3 < 30, < A, < 0.8, < mq, A^, x, niA < 1000 GeV, 
and Q < At < 2000 GeV, using a random-number generating function. The number of 
points we explore in this region of the parameter space is 2 x 10^. Among those 2 x 10^ 
points randomly distributed in the parameter space, only 14571 points are physically 
allowed. At these 14571 points of the parameter space, we consider the neutral Higgs 
productions at the ILC, the future high-energy e^e~ colliding machine, with the proposed 
center of mass energy of 500 GeV. 

A. HIGGS MASS 

The mass of the lightest neutral Higgs boson, rrih^, is calculated at each of the physically 
allowed points we explore in the prescribed region, and the masses of heavier neutral 
Higgs bosons as well. The result is shown in Fig. 1, where rrth^ is plotted against tan/3. 
One can find that the upper bound on the lightest neutral Higgs boson mass increases 
as tan (3 increases, and it is nearly saturated at about 150 GeV as tan (3 approaches to 
30. The behavior of the upper bound on the lightest neutral Higgs boson mass in the 
NMSSM with explicit CP violation is somewhat similar to that of the MSSM with no CP 
violation. It may be contrasted against the case of no CP violation in the Higgs sector of 
the NMSSM, where the upper bound on the lightest neutral Higgs boson mass becomes 
maximal at around tan/? ~ 2.5 and decreases slightly as tan/3 further increases, for a 
wide region in the parameter space, mainly because of the Higgs singlet effects. 

B. HIGGS PRODUCTION 

The most important channels for the production of neutral Higgs bosons in high energy 
e"'"e~ collisions are: the Higgs-strahlung process e+e" Zhi, the WW fusion process 
e"'"e~ — > VfAJfhi, and the ZZ fusion process e^e~ — e^e~hi [i = 1-5). We denote the 
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production cross sections for each of these processes as af{mh^), aj^{mh-), and afirrih-), 
respectively. These production cross sections in each channel are related to the production 
cross section for the SM Higgs boson in the corresponding channel as [20] 

where Ri {i = 1-5) are defined as 

Ri = cos fSOu + sin f302i , (15) 

with Oij being the elements of the orthogonal matrix that diagonalizes the 5x5 mass 
matrix of the neutral Higgs bosons. They satisfy < i?f < 1 (i = 1-5) and the orthogonal 
condition of I]f=i -R? = 1- Note that they represent collectively all relevant parameters of 
the NMSSM at the one-loop level, including the CP violating phases. 

Let us first consider the Higgs-strahlung process for the neutral Higgs productions. 
At each physically allowed point, we calculate the mass of the participating neutral Higgs 
boson, and its production cross section via the Higgs-strahlung process. We repeat this 
job for all five neutral Higgs bosons, to obtain af{mhi) {i = 1-5) at each of 14571 points. 
Now, we denote by ao{hZZ) the largest production cross section among them, namely, 

aoihZZ) = max[o-f , af , , , a^] . (16) 

Thus, at every physically allowed point of the parameter space of the NMSSM, ao{hZZ) 
is obtained. The meaning of ao{hZZ) is that, at a given point of the parameter space, 
at least one of the five neutral Higgs bosons in the NMSSM may be produced via the 
Higgs-strahlung process with that cross section. 

In Fig 2a, we plot ao{hZZ) against m/^ for the future high energy e"'"e~ collider with 
y/s = 500 GeV. The dashed curve in the figure represents the production cross section 
of the SM Higgs boson via the Higgs-strahlung process [20,21]. Among ao{hZZ), we 
may select the smallest one and denote it as . The significance of is quite clear; 
at least one of the five neutral Higgs bosons in the NMSSM may be produced via the 
Higgs-strahlung process with a cross section larger than cx^, regardless of the parameter 
values. In this sense, it may be regarded as the absolute minimum of the production cross 
section of a neutral Higgs boson in the NMSSM via the Higgs-strahlung process, with 
explicit CP violation at the one-loop level. From Fig. 2a, one can see that ~ 12 fb. 

Next, we consider the neutral Higgs productions via the WW fusion process and the 
ZZ fusion process. For the WW fusion process, we determine (Jo(/iz/z/) and , which are 
respectively analogous to ao{hZZ) and for the Higgs-strahlung process. We show the 
result for aolhuu) in Fig. 2b, plotted against rrih^, and we obtain ~ 15 fb. Also, for 
the ZZ fusion process, we determine ao{he~^e~) and , which are respectively analogous 
to a^iJiZZ) and for the Higgs-strahlung process. We show the result for (To(/ie+e~) 
in Fig. 2c, plotted against m/i^, and we obtain ~ 1.5 fb. The dashed curves in Figs. 
2b and 2c represent respectively the production cross sections of the SM Higgs boson via 
the WW fusion and the ZZ fusion processes. 

We find that, interestingly, all of ~ 12 fb, ~ 15 fb, and ~ 1.5 fb occur at an 
identical point of the parameter spce, where the values of the parameters are ^ </)f ^ 7r/6, 
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tan/3 = 20, mq = 766 GeV, A = 0.088, k = 0.119, Ak = 86 GeV, x = 699 GeV, At = 1350 
GeV, and rriA = 139 GeV. We may regard this point of the parameter space with these 
values as the most difficult point for the neutral Higgs productions in e~^e~ collisions with 
^/s = 500 GeV. We would like to investigate the region around this point in more detail. 

Let us study a narrow region defined as < 0j = < vr, < tan/3 < 20, mg = 766 
GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, A = 1350 GeV, and = 
139 GeV. It is in fact a two-dimensional subspace in the (0, tan/3)-plane. We calculate 
ao{hZZ) in this restricted region. The result is shown in Fig. 3a. Likewise, the result 
for ao^huu) and ao{he^e~) are respectively shown in Figs. 3b and 3c. One can see that 
indeed ~ 12 fb, ~ 15 fb, and ~ 1.5 fb at the point (0, tan/5) = (7r/6,20) in 
the (0, tan/?)-plane, indicated by the arrow, in Figs. 3a, 3b, and 3c, respectively, for the 
neutral Higgs productions in e~^e~ collisions with = 500 GeV. 

Now, we study the effect of the CP phases on the neutral Higgs productions in e+e^ 
collisions with ^/s = 500 GeV. The parameters are set as < </)t = </> < vr, tan/3 = 20, 
mq = 766 GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, At = 1350 GeV, and 
rriA = 139 GeV. Thus, we are investigating a one-dimensional subspace of the line in 
direction. 

In this one-dimensional subspace of the parameter space, we calculate the masses of 
the five neutral Higgs bosons in the present model, in order to see the dependencies of the 
neutral Higgs masses upon the CP phase. The results are shown in Fig. 4a, where the 
five neutral Higgs boson masses are plotted as functions of 0. One can notice in Fig. 4a 
that the five neutral Higgs boson masses may vary up to about 5 GeV, as the CP phase 
changes within the given range. 

The coefficients Ri [i = 1-5), which relate the neutral Higgs production cross sections 
in the present model with the production cross sections for the SM Higgs boson, also 
depend on the parameters of the present model. We plot in Fig. 4b Ri as functions of 
0. All Ri fluctuate quite significantly as changes within the given range. These wild 
behaviors of Ri indicate that the production cross sections of the neutral Higgs bosons 
would also fluctuate largely against the CP phase. The production cross sections of the 
five neutral Higgs bosons are displayed as functions of 0, via the Higgs-strahlung, WW 
fusion, and ZZ fusion processes, respectively, in Figs. 4c, 4d, and 4e. As might be 
expected, all of the neutral Higgs production cross sections fluctuate largely against 0. 
In fact, the patterns of fluctuation in the neutral Higgs productions is very similar to 
the fluctuation patterns in Ri. This is mainly because the dependencies of the masses of 
the neutral Higgs bosons upon the CP phase are relatively negligible as compared to the 
dependencies of Ri upon the CP phase. 

The three figures. Figs. 4c, 4d, and 4e, are useful to determine a^ifiZZ), aolhuu), and 
ao{he~^e~), as well as alf , , and , respectively. Consider for example the neutral 
Higgs productions via the Higgs-strahlung process in Fig. 4c, where all of , , (t^ ■, 

, and are shown. Since aoihZZ) is defined as ao{hZZ) = msix[a^ , (t|^, , af, a^], 
we may identify the envelope of the five curves in Fig. 4c as ao{hZZ). Thus, in Fig. 
4c, the curve for ao{hZZ) is somewhat W-shaped. One may notice that the minimum 
of ao{hZZ), which is defined as cx^, occurs at ~ G/n with cr^ ~ 12 fb. Likewise, one 
may also notice in Figs. 4d and 4e respectively that ~ 15 fb for the neutral Higgs 
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productions via the WW fusion process and o"^ ~ 1.5 fb for the neutral Higgs productions 
via the ZZ fusion process. The results obtained from Figs. 4c, 4d, and 4e are consistent 
with the results from Figs. 2a, 2b, and 2c, respectively. 



C. HIGGS DECAY 

Now, let us study the decay modes of the five neutral Higgs bosons hj (j = 1-5) in the 
NMSSM with explicit CP violation. The partial decay width of a neutral Higgs boson hj 
into a pair of down or up quarks is given as 



dd) 



uu] 



327rmfj/ 



rd 



327rm2 



To, 



W 



cos^ P 
sin^ P 



rd 



tan^ (30 



i3 



Tu) + coe (30% 



(17) 



where = 3 is the color factor of quarks, and Tq = Arii?^/m\, [q = d,u). The color factor 
of leptons is 1. The partial decay width of a neutral Higgs boson hj into a pair of charged 
leptons is the same as T{hj — > dd) with Cf = 1 for the color factor of leptons. 

The partial decay width of hj into a pair of gauge bosons may be obtained from that 
of the SM Higgs boson, through a relation given by 



Tih, VV) 



Rj ^suihj 



VV) 



(18) 



where rsM(^i — > VV) is the decay width of the SM Higgs boson into a pair of gauge 
bosons. The produced gauge bosons may be either real or virtual. 
The partial decay width of hj into a pair of photons is given as 



T{hj 77) 



am 



hi 
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+ 



a mi 
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(19) 



where and ei, are respectively the electric charges of top and bottom quarks, tw = 
im^y/ml^., and Sq{Tq), Pq{Tq) (g = t,6), and Sw{tw) are respectively the form factors of 
the scalar, pseudoscalar, and gauge bosons, given as [10,22] 



Sq^Tq) 



Irqil 
3 



^)f{rq)) 



[2 + 3tw + 3rH/(2 



where the function / is defined as 
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(20) 
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Lastly, the partial decay width of hj into a pair of gluons is given as 




71 




' 2887r3m^ |tan/30,3n(r6) + cot/30,3Pt(n)|' 

J- + -fiVj^ . 

71 



(22) 



where as{mh ) is the coupling coefficient of the strong interactions, Np is the number of 
the quark flavors lighter than hj, and 




Now, we calculate the total decay width of the neutral Higgs bosons in the NMSSM 
with explicit CP violation, in the parameter range set for Fig. 4a, namely, < (/> < tt 
{(j) = (j)^ = </,f,), tan/? = 20, tuq = 766 GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 
GeV, At = 1350 GeV, and rriA = 139 GeV. Let us remind the reader that the neutral 
Higgs bosons are most difficult to be discovered in this parameter range. The masses of all 
neutral Higgs bosons lie within the range of 120 to 160 GeV, as shown in Fig. 4a. Thus, 
those decay modes listed above are the dominant channels for the Higgs boson research. 
We assume that the total decay widths of hj is given by 



r{hj) = r{hj -> bb) + r{hj T+T-) + r{hj cc) + r{hj ss) 

+ T{hj ^ W+W-) + r(/i, ^ ZZ) + T{hj ^ 77) + T{hj ^ gg) . (24) 



In Fig. 5a, the total decay widths of the five neutral Higgs bosons are plotted as functions 
of the phase 0. Note that the total decay widths fluctuate widely against 0. The results 
shown in Fig. 5a may be compared to the case of the SM, where the total decay width 
of the SM Higgs boson is obtained between 0.003535 GeV and 0.08613 GeV as the SM 
Higgs boson mass increases from 120 GeV to 160 GeV. On the other hand, the total decay 
width of hi in the present model may be as large as about 0.63 GeV, according to Fig. 
5a. 

In Figs. 5b, 5c, 5d, 5e, and 5f, the branching ratios of each of the flve neutral Higgs 
bosons are shown as functions of 0, in the same parameter range set for Fig. 4a or Fig. 
5a. One may easily observe, in these flgures, that for the whole range of < < vr, 
T{hj bb) is dominant over other partial decay widths, for all flve neutral Higgs bosons, 
except for Fig. 5e, where the branching ratio of h^ into W^W~ is larger than that of h^ 
into hb, in a narrow range near ~ 27r/3. The decay widths of the flve neutral Higgs 
bosons into r"*"r~ are also important, except /14, in the present model, at least for larger 
tan (5. Note that we set tan /3 = 20 in our numerical calculations. 
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One may also observe that for the most part of the range the branching ratio of 
hj into ss is larger than that of into cc, for all five neutral Higgs bosons except /14. We 
have a large branching ratio of hj into ss, mainly because we set tan/3 = 20, since the 
decay width of a neutral Higgs boson into a pair of down-type quarks, as well as that into 
a pair of charged leptons, increases as tan f3 increases. This observation is interesting as 
compared to the case of the SM, where the branching ratio of the SM Higgs boson into 
cc is larger than that into ss because the coupling strength of the SM Higgs boson with 
a fermion is proportional to the mass of the fermion. 

IV. CONCLUSIONS 

We have studied the Higgs sector of the NMSSM with explicit CP violation, at the one- 
loop level where the radiative corrections due to the quarks and squarks of the third 
generation are taken into account. The CP-odd tadpole minimum conditions at the one- 
loop level are analytically obtained, and the mass matrix for the five neutral Higgs bosons 
are calculated. Unlike the CP conserving case in the Higgs sector of the NMSSM, the 
upper bound on the mass of the lightest neutral Higgs boson in the present model is found 
to increase as tan f3 increases, in the parameter region set as < 0, 0^ < tt, < A, < 0.8, 
< mg, Ak, X, TUA < 1000 GeV, and < At < 2000 GeV. 

The production cross sections of the five neutral Higgs bosons, via the Higgs-strahlung 
process, the WW fusion process, and the ZZ fusion process, are calculated in the pa- 
rameter region set as < 0, 0( < vr, < A, /c < 0.8, < niQ, Ak,x,mA_ < 1000 GeV, 
and < At < 2000 GeV. It is found that, in this parameter range, at least one of the 
five neutral Higgs bosons may be produced at the ILC with ^/s = 500 GeV with cross 
sections larger than 12, 15, and 1.5 fb, respectively, via the Higgs-strahlung process, the 
WW fusion process, and the ZZ fusion process. These numbers suggest that it is possible 
for the ILC with ^/s = 500 GeV to test the present model. 

Further, the effects of the CP phase on the productions of the five neutral Higgs bosons 
via the Higgs-strahlung process, the WW fusion process, and the ZZ fusion process, are 
studied for the parameter region set as < < vr (0 = 0t = 0;,), tan/3 = 20, rriQ = 766 
GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, At = 1350 GeV, and tua = 139 
GeV. This is a one-dimensional subspace in direction of the parameter space where it is 
found that the production cross sections for at least one of the five neutral Higgs bosons 
are minimum. In this one- dimensional parameter space of < < vr, the cross sections 
of the five neutral Higgs bosons are found to be seriously dependent on 0, regardless of 
via which process they are produced. Their large fluctuations against suggest that the 
explicit CP violation in the NMSSM through is indeed very crucial for the study of the 
Higgs sector in the present model. 

The decay widths of the five neutral Higgs bosons into various particles are also studied 
in the parameter region set as < < tt (0 = 0t = 0^), tan/3 = 20, mg = 766 GeV, 
A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, At = 1350 GeV, and niA = 139 
GeV. It is found that the present model predicts somewhat different patterns of neutral 
Higgs boson decays from those of the SM or those of the NMSSM with no CP violation. 
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Whereas the branching ratio of the SM Higgs boson into cc is larger than that into ss, 
the branching ratio of hj in the present into ss is larger than that into cc, for the most 
part of the range, for all five neutral Higgs bosons except h^. 

In conclusion, the scenario of explicit CP violation in the Higgs sector of the NMSSM 
exhibits interesting behaviors of the present model and thus is worthwhile for further 
study. 
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FIGURE CAPTIONS 



Fig. 1: The mass of the hghtest neutral Higgs boson is plotted against tan/?, for < 
0, 04 < TT, < A, A; < 0.8, < mg, A^, x, tua < 1000 GeV, and < At < 2000 GeV. 

Fig. 2(a): The plot of ao{hZZ) against rrih^ in e~^e~ collisions with ^/s = 500 GeV, for 
< (j),(f)t < rr, 2 < tan/5 < 30, < A, A; < 0.8, < tuq, Ak,x,mA < 1000 GeV, and 
< At < 2000 GeV. The dashed curve represents the production cross section of the SM 
Higgs boson via the Higgs-strahlung process. The minimum value of aQ{hZZ) is about 
12 fb, represented by a point in the lower right corner at about m/j^ ~ 127 GeV. 

Fig. 2(b): The plot of aoi^huu) against m/^ in e+e" collisions with y/s = 500 GeV, for 
the same parameter ranges as Fig. 2a. The dashed curve represents the production cross 
section of the SM Higgs boson via the WW fusion process. The minimum value of ao{hujj) 
is about 15 fb, represented by a point in the lower right corner at about nih^ ~ 127 GeV. 

Fig. 2(c): The plot of ao{he~^e~) against m^^ in e+e^ collisions with y/s = 500 GeV, 
for the same parameter ranges as Fig. 2a. The dashed curve represents the production 
cross section of the SM Higgs boson via the ZZ fusion process. The minimum value 
of (XQ{he~^e~) is about 1.5 fb, represented by a point in the lower right corner at about 
ruh, ~ 127 GeV. 

Fig. 3(a): The lego plot of ao{hZZ) in the (tan /3, 0)-plane for the Higgs production via 
the Higgs-strahlung process in e~^e~ collisions with ^/s = 500 GeV, for (pt = </>, rnq = 766 
GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, A = 1350 GeV, and = 139 
GeV. The minimum value of aQ{hZZ) is about 12 fb, which occurs at a dip in the forefront 
of the (tan /3, 0)-plane at (tan/3, 0) = (20,3/7r), indicated by an arrow. 

Fig. 3(b): The lego plot of aQihi/v) in the (tan /?, 0)-plane for the Higgs production via 
the WW fusion process in e"'"e~ collisions with a/s = 500 GeV, for the same parameters 
as Fig. 3a. The minimum value of (To(/ii^z/) is about 15 fb, which occurs at the same point 
of the (tan /?, {/))-plane as Fig. 3a. 

Fig. 3(c): The lego plot of ao{he^e~) in the (tan/?, </))-plane for the Higgs production via 
the ZZ fusion process in e^e~ collisions with \/s = 500 GeV, for the same parameters 
as Fig. 3a. The minimum value of ao{he^e~) is about 1.5 fb, which occurs at the same 
point of the (tan /5, (/))-plane as Fig. 3a. 

Fig. 4(a): The masses of the five neutral Higgs bosons are plotted as functions of 0, 
where the parameters are set as (f>t = 0, tan/? = 20, mg = 766 GeV, A = 0.088, k = 0.12, 
Ak = 86 GeV, x = 699 GeV, At = 1350 GeV, and = 139 GeV. 
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Fig. 4(b): Ri {i = 1-5) as functions of 0, for the same parameters as Fig. 4a. 

Fig. 4(c): The production cross sections of hi {i = 1-5) via the Higgs-strahlung process 
in e+e" colhsions with ^/s = 500 GeV, as functions of </>, for the same parameters as Fig. 
4a. 

Fig. 4(d): The production cross sections of hi {i = 1-5) via the WW fusion process in 
e~^e~ colhsions with y/s = 500 GeV, as functions of 0, for the same parameters as Fig. 
4a. 

Fig. 4(e): The production cross sections of hi {i = 1-5) via the ZZ fusion process in e+e~ 
coUisions with ^/s = 500 GeV, as functions of (p, for the same parameters as Fig. 4a. 

Fig. 5(a): The total decay widths of the five neutral Higgs bosons as functions of 0, for 
the same parameters as Fig. 4a. 

Fig. 5(b): The branching ratios of hi as functions of 0, for the same parameters as Fig. 
4a. 

Fig. 5(c): The branching ratios of /i2 as functions of (p, for the same parameters as Fig. 
4a. 

Fig. 5(d): The branching ratios of h^ as functions of 0, for the same parameters as Fig. 
4a. 

Fig. 5(e): The branching ratios of h^ as functions of 0, for the same parameters as Fig. 
4a. 

Fig. 5(f): The branching ratios of h^ as functions of 0, for the same parameters as Fig. 
4a. 
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Figure 1: The mass of the hghtest neutral Higgs boson is plotted against tan/9, for 
< 0, 04 < TT, < A, A; < 0.8, < tuq, Ak, x, tua < 1000 GeV, and < A < 2000 GeV. 
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Figure 2a: The plot of ao{hZZ) against m^^ in e"*"e~ collisions with ^/s = 500 GeV, for 
< 0, 0t < TT, 2 < tan/5 < 30, < A, < 0.8, < mq, Ak,x,mA < 1000 GeV, and 
< At < 2000 GeV. The dashed curve represents the production cross section of the SM 
Higgs boson via the Higgs-strahlung process. The minimum value of ao{hZZ) is about 
12 fb, represented by a point in the lower right corner at about rrih^ ~ 127 GeV. 
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Figure 2b: The plot of aoi^huu) against rrih^ in e^e^ collisions with ^/s = 500 GeV, for 
the same parameter ranges as Fig. 2a. The dashed curve represents the production cross 
section of the SM Higgs boson via the WW fusion process. The minimum value of aoihuu) 
is about 15 fb, represented by a point in the lower right corner at about rrih^ ~ 127 GeV. 
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Figure 2c: The plot of ao{he~^e~) against nih^ in e~^e~ collisions with y/s = 500 GeV, 
for the same parameter ranges as Fig. 2a. The dashed curve represents the production 
cross section of the SM Higgs boson via the ZZ fusion process. The minimum value 
of ao{he~^e~) is about 1.5 fb, represented by a point in the lower right corner at about 
nih, ~ 127 GeV. 
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Figure 3a: The lego plot of ao{hZZ) in the (tan /3, 0)-plane for the Higgs production via 
the Higgs-strahlung process in e~^e~ collisions with ^/s = 500 GeV, for (pt = <f), friq = 766 
GeV, A = 0.088, k = 0.12, Ak = 86 GeV, x = 699 GeV, At = 1350 GeV, and tua = 139 
GeV. The minimum value of ao{hZZ) is about 12 fb, which occurs at a dip in the forefront 
of the (tan /3, 0)-plane at (tan/3, 0) = (20,3/7r), indicated by an arrow. 
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Figure 3b: The lego plot of aoihuu) in the (tan /3, 0)-plane for the Higgs production via 
the WW fusion process in e~^e~ collisions with ^/s = 500 GeV, for the same parameters 
as Fig. 3a. The minimum value of aolhuu) is about 15 fb, which occurs at the same point 
of the (tan (3, 0)-plane as Fig. 3a. 
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Figure 3c: The lego plot of ao{he~^e~') in the (tan/3, 0)-plane for the Higgs production via 
the ZZ fusion process in e^e~ collisions with ^/s = 500 GeV, for the same parameters 
as Fig. 3a. The minimum value of cro(/ie~^e~) is about 1.5 fb, which occurs at the same 
point of the (tan /3, 0)-plane as Fig. 3a. 
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Figure 4a: The masses of the five neutral Higgs bosons are plotted as functions of 
where the parameters are set as (pt = (p, tan/5 = 20, mq = 766 GeV, A = 0.088, k = 0.12, 
Ak = 86 GeV, x = 699 GeV, A = 1350 GeV, and tua = 139 GeV. 
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Figure 4c: The production cross sections of hi {i = 1-5) via the Higgs-strahlung process 
in e~^e~ colhsions with ^/s = 500 GeV, as functions of 0, for the same parameters as Fig. 
4a. 
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Figure 4d: The production cross sections of hi {i = 1-5) via the WW fusion process in 
e'^e~ colhsions with ^/s = 500 GeV, as functions of 0, for the same parameters as Fig. 
4a. 
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Figure 5b: The branching ratios of hi as functions of 0, for the same parameters as Fij 
4a. 
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Figure 5f: The branching ratios of as functions of 0, for the same parameters as Fig. 
4a. 
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